ID  6 
ID  g 

m~m  CO 

<s 


DAVID  W.  TAYLOR  NAVAL  SHIP 
RESEARCH  AND  DEVELOPMENT  CENTER 

Bethesda,  Maryland  20084-5000 


MEASURED  SIMILARITY  PROPERTIES  OF  EDDY  VISCOSITY 
AND  MIXING  LENGTH  IN  THREE-DIMENSIONAL 
TURBULENT  STERN  FLOWS 


Thomas  T.  Huang 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  IS  UNLIMITED. 


Presented  at  the 

Second  International  Symposium  on  Ship 
Viscous  Resistance  1985,  SSPA 
Goteborg,  Sweden 
18-20  March  1985 


SHIP  PERFORMANCE  DEPARTMENT 
RESEARCH  AND  DEVELOPMENT  REPORT 


electei 

JUL1  719851 


May  1985 


DTNSRDC -85/038 


85  7  03  006 


MAJOR  DTNSRDC  ORGANIZATIONAL  COMPONENTS 


uh'0  «<l7-440 


NDW  DTNSRDC  39«0/43h  Z'SCM 


UNCLASSIFIED 


SECURITY  classification  of  this  page 


lj  REPORT  SECURITY  CLASSIFICATION 
UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
DTNSRDC-85/038 


3  DISTRIBUTION/ AVAILABILITY  OF  REPORT 
APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  IS  UNLIMITED. 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


ba  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
David  W.  Taylor  Naval  Ship  {if  applicable) 

R  and  D  Center  1542 


6c  ADDRESS  (C/fy,  State  and  HP  Code) 

Bethesda,  Maryland  20084-5000 


7b  ADDRESS  (City.  State,  and  ZIP  Code) 


Ra  NAME  OF  FUNDING  <  SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL 
(If  applicable) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


<  ADDRESS  (City.  State  and  ZIP  Code) 

'0  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 

PROJECT 

TASK 

ELEMENT  NO 

NO 

NO 

61152N 

ZR0230101 

WORK  UNIT 
ACCESSION  NO 


’l  NTiE  (Include  Security  Classification) 

MEASURED  SIMILARITY  PROPERTIES  OF  EDDY  VISCOSITY  AND  MIXING  LENGTH  IN  THREE-DIMENSIONAL 
TURBULENT  STERN  FLOWS 


'2  -ORSON A.  AuThOR(S) 
Thomas  T.  Huang 


•  >a  Tyre  OF  REPORT 
Finn  L 


1  3b  TIME  COVERED 
FROM  TO 


14  DATE  OF  REPORT  (Year,  Month.  Day)  15  PAGE  COUNT 

1985  May  31 


•b  v  PR-  f  VENTARY  NOTATION 

Presented  at  the  Second  International  Symposium  on  Ship  Viscous  Resistance  1985,  SSPA 
Coteborg,  Sweden,  18-20  March  1985 


'8  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 
Thick  turbulent  stern  flows,  turbulence  modeling, 
eddy  viscosity,  mixing  length.^ 


'9  a-IS'Rai  t  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

"* A  comprehensive  set  of  experimental  turbulence  data  is  presented  for  the  stern  region  of 
simple  three-dimensional  models  having  1:1,  2:1,  and  3:1  elliptic  transverse  cross 
sections.  Direct  measurements  of  Reynolds  stresses,  turbulence  kinetic  energy,  and 
mean  velocity  profiles  are  made  by  "X"  hot-film  anemometry.  The  values  of  eddy  viscosity 
and  mixing  length  are  obtained  from  the  measured  Reynolds  stresses  and  the  mean  velocity 
gradient.  The  experimentally-determined  distributions  of  eddy  viscosity  and  mixing 
length  in  the  thick  stern  boundary  layers  are  found  to  be  smaller  than  values  which  have 
been  proposed  for  thin  boundary  layers.  Because  eddy  viscosity  and  mixing  length  models 
play  an  important  role  in  stern  boundary-layer  calculations,  a  revised  empirical  mixing 
length  model,  which  is  assumed  to  be  proportional  to  the  square-root  of  the  effective 
turbulence  area,  is  proposed.  The  similarity  properties  of  eddy  viscosity  and  mixing 


C  o.  w 


.’0  >Sr"'9-,rON  AVAILABILITY  OF  ABSTRACT 

.NCt  ASSlFiED  UNLIMITED  □  SAME  AS  RPT  Q  DTlC  USERS 


'.M  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Thomas  T.  Huan 


(Continued  on  reverse  side) 


21  ABSTRACT  SECURITY  CLASSIFICATION 
UNCLASSIFIED 


22b  TELEPHONE  (Include  Area  Code)  22 c  OFFICE  SYMBOL 

(202)  227-1325  Code  1542 


DO  FORM  1473,  84  MAR  83  APR  edition  may  be  used  until  exhausted 

All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE_ 

UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


l.t  RFPOH1  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


Zb  DECLASSIFICATION/  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

DTNSRDC-85/038 


6a  NAME  OF  PERFORMING  ORGANIZATION 
David  W.  Taylor  Naval  Ship 
R  and  D  Center 


6c  ADDRESS  (City,  Sfdfe.  and  /IPCode) 

Bethesda,  Maryland  20084-5000 


6b  OFFICE  SYMBOL 
(If  applicable) 

1542 


Ra  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 


■\DDRFSS  (Cify,  Sfdfe,  and  ZIP  Code) 


8b  OFFICE  SYMBOL 
(If  applicable)  ■ 


pb  RESTRICTIVE  MARKINGS 


1  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  IS  UNLIMITED. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7a  NAME  OF  MONITORING  ORGANIZATION 


7b  ADDRESS  (C/ty.  State,  and  ZIP  Code) 


PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


I  WORK  UNIT 
ACCESSION  NO 


61152N  ZR0230101 


TIDE  ( include  Security  Classification) 

MEASURED  SIMILARITY  PROPERTIES  OF  EDDY  VISCOSITY  AND  MIXING  LENGTH  IN  THREE-DIMENSIONAL 
TURBULENT  STERN  FLOWS 


1Z  PERSONAL  AUTHOR(S) 

Thomas  T.  Huang 


1  Id  TYPE  OF  REPORT 
Final 


13b  TIME  COVERED 
FROM  TO 


14  DATE  OF  REPORT  (Year,  Month.  Day)  15  PAGE  COUNT 

1985  May  31 


"j  SUPPLEMENTARY  NOTATION 

Presented  at  the  Second  International  Symposium  on  Ship  Viscous  Resistance  1985,  SSPA 
Gotcborg,  Sweden,  18-20  March  1985 


18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number)  \ 
'^Thick  turbulent  stern  flows,  turbulence  modeling, 
eddy  viscosity,  mixing  length,^- 


COSATI  CODES _ 

GROUP  I  SUBGROUP 


ri» 


abstract  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  comprehensive  set  of  experimental  turbulence  data  is  presented  for  the  stern  region  of 
simple  three-dimensional  models  having  1:1,  2:1,  and  3:1  elliptic  transverse  cross 
sections.  Direct  measurements  of  Reynolds  stresses,  turbulence  kinetic  energy,  and 
mean  velocity  profiles  are  made  by  "X"  hot-film  anemometry.  The  values  of  eddy  viscosity 
and  mixing  length  are  obtained  from  the  measured  Reynolds  stresses  and  the  mean  velocity 
grad  Lent.  The  experimentally-determined  distributions  of  eddy  viscosity  and  mixing 
length  Ln  the  thick  stern  boundary  layers  are  found  to  be  smaller  than  values  which  have 
been  proposed  for  thin  boundary  layers.  Because  eddy  viscosity  and  mixing  length  models 
p  Lay  an  important  role  in  stern  boundary-layer  calculations,  a  revised  empirical  mixing 
length  model,  which  Ls  assumed  to  be  proportional  to  the  square-root  of  the  effective 
turbulence  area,  is  proposed.  The  similarity  properties  of  eddy  viscosity  and  mixing 

~'a  Ct  f<A  M  LU- I  ,  i  . 

(,  _ (Continued  on  reverse  side) 


70  DISTRIBUTION  /  AVAILABILITY  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

fcluNCLASSIFIEDfUNLIMITED  □  SAME  AS  RPT  □  DTIC  USERS  UNCLASSIFIED _ 


.Vj  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 

Dr.  Thomas  T.  Huang  (202)  227-1325 _ Code  1542 _ 


DO  FORM  1 473,  84  MAH  83  APR  edition  may  be  med  until  e»hau»ted  _ SECURITY  CLASSIFICATION  OF  rm$  PA£f _ 

All  other  edition*  are  obtolete  UNCLASSIFIED  ~  ' 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


(Block  19  Continued) 

Length  can  be  incorporated  into  the  three-dimensional  boundary  layer  computation 
methods  using  simple  turbulence  modeling  or  can  be  used  as  an  initial  input  to  the 
two-equation  (k-e)  turbulence  modeling  method.  The  measured  turbulence  structural 

parameters  -u^vV(2K)  and  v  t/(/k£)  are  found  to  be  almost  constant  across  the  entire 

thick  stern  boundary  layer.  The  measured  proportionality  constants  are  in  agreement 
with  those  derived  from  the  near-wall  turbulence.  The  distributions  of  turbulence 
kinetic  energy  k  and  the  turbulence  energy  dissipation  £  can  be  approximately 
estimated  once  the  distributions  of  mixing  length  ?,  and  hence  the  eddy  viscosity 

v  =  Du  /Dn  and  the  dominant  Reynolds  stress  -u^v 7  =  v  Du  /Dn  are  known,  where  Du  /Dn 
t  x  t  x  x 

is  the  primary  velocity  gradient  of  the  flow. 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


TABLE  OF  CONTENTS 


LIST  OF  FIGURES  . 

TABLE  . 

ADMINISTRATIVE  INFORMATION  . 

ABSTRACT  . 

NOTATION  . 

1.  INTRODUCTION  . 

2.  MODELS  AND  EXPERIMENTAL  TECHNIQUES  . 

3.  MEASURED  AND  DERIVED  TURBULENT  CHARACTERISTICS 

OF  THICK  STERN  FLOWS  . 

3.1  Measured  Reynolds  Stresses  . 

3.2  Measured  Turbulence  Kinetic  Energy  k  . 

3.3  Derived  Eddy  Viscosity,  Mixing  Length,  and 
Turbulence  Kinetic  Energy  Dissipation  .  .  .  . 

4.  COMPARISON  OF  SIMPLE  THIN  BOUNDARY  LAYER  TURBULENCE 

MODELS  WITH  THE  MEASURED  DATA  . 

4.1  Comparison  of  the  Bradshaw  et  al.  (7)  Mixing 

Length  Model  with  Present  Data  . 

4.2  Comparison  of  the  Cebeci  and  Smith  (8)  Eddy 

Viscosity  Model  with  Present  Data  . 

5.  "ZONAL"  TURBULENT  SIMILARITY  PROPERTIES  OF  THICK 

STERN  FLOWS . . 

5.1  Similarity  Property  of  Mixing  Length  for 

Turbulent  Stern  Flows  . 

5.2  Measured  Turbulence  Structural  Parameters  .  . 

6.  CONCLUSION  . 

ACKNOWLEDGMENTS  . 

REFERENCES  . 


LIST  OF  FIGURES 


Page 


1  -  Schematic  of  the  Three  Sterns  Showing  Measurement  Stations  .  14 

2  -  Two  Elliptical  Models  Mounted  in  a  Wind  Tunnel .  14 

3  -  Measured  Distributions  of  Reynolds  Stresses  (Model  1:1)  .  15 

4  -  Measured  Turbulence  Kinetic  Energy  .  16 

5  -  Measured  Distribution  of  Mixing  Length  Parameter  .  17 

6  -  Measured  Distribution  of  Eddy  Viscosity  .  18 

7  -  Representing  the  Mixing  Length  of  the  Stern  Flow  by  the 

Square  Root  of  Turbulence  Area .  19 

8  -  Similarity  Property  of  Mixing  Length  for  Turbulent  Stern  Flows  .  20 

9  -  Schematic  of  Mixing  Length  across  the  Boundary  Layer  .  .  21 

10  -  Measured  Distribution  of  Turbulence  Structure  Parameter  .  21 

11  -  The  Constancy  of  the  Measured  Values  of  V  //k  &  across  the 

Thick  Boundary  Layers .  22 


Table  1  -  Offset  of  DTNSRDC  Axisymmetric  Model  1 .  23 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  paper  was  funded  under  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center's  Independent  Research  Program,  Program  Element 
61152N,  Task  Area  ZR0230101,  and  Work  Unit  1542-103. 


\  \  ^  - — ■ — *: — ■v  -•  v*-  -.  ^  -  V’-.  — ■ 


MEASURED  SIMILARITY  PROPERTIES  OF  EDDY  VISCOSITY 
AND  MIXING  LENGTH  IN  THREE-DIMENSIONAL  TURBULENT 
STERN  FLOWS 

T.  T.  HUANG 

David  W.  Taylor  Naval  Ship  Research  &  Development  Center 
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ABSTRACT 

A  comprehensive  set  of  experlmentsl  turbulence  date  Is  presented  for  the  stern  region  of 
simple  three-dimensional  models  having  1:1,  2:1,  and  3:1  elliptic  transverse  cross  sections. 
Direct  measurements  of  Reynolds  stresses,  turbulence  kinetic  energy,  and  mean  velocity  pro¬ 
files  are  made  by  "X"  hot-fllm  anemometry.  The  values  of  eddy  viscosity  and  mixing  length 
are  obtained  from  the  measured  Reynolds  stresses  and  the  mean  velocity  gradient.  The 
experimentally-determined  distributions  of  edd,  viscosity  and  mixing  length  in  the  thick 
stern  boundary  layers  are  found  to  be  smaller  than  values  which  have  been  proposed  for  thin 
boundary  layers.  Because  eddy  viscosity  and  mixing  length  models  play  an  important  role  in 
stern  boundary- layer  calculations,  a  revised  empirical  mixing  length  model,  which  la  assumed 
to  be  proportional  to  the  square-root  of  the  effective  turbulence  area.  Is  proposed.  The 
similarity  properties  of  eddy  viscosity  and  mixing  length  can  be  Incorporated  Into  the 
three-dimensional  boundary  layer  computation  methods  using  simple  turbulence  modeling  or  can 
be  used  as  an  Initial  Input  to  the  two-equation  (k-e)  turbulence  modeling  method.  The 

measured  turbulence  structural  parameters-u  v  /(2K)  and  vt/(/kt)  are  found  to  be  almost 
constant  across  the  entire  thick  stern  boundary  layer.  The  measured  proportionality 
constants  are  In  agreement  with  those  derived  from  the  near-wall  turbulence.  The  distribu¬ 
tions  of  turbulence  kinetic  energy  k  and  turbulence  energy  dissipation  E  can  be  approxima¬ 
tely  estimated  once  the  distributions  of  mixing  length  i  and  hence  the  eddy  viscosity 

Vt  -  3ux/3n  and  the  dominant  Reynolds  stress  -u  v’  -  vt  3ux/3n  are  known,  where  3ux/3n  is 
the  primary  velocity  gradient  of  the  flow. 

NOTATION 

A  Van  Drlest's  damping  factor  (see  Equation  6) 

A(x)  Cross  turbulence  area,  A(x)/x  *  [  (a+0.  65a)(b+0.  66{,)  -  (a+ca)(b+Eb)  I 
a,b  Length  of  the  malor  and  minor  elliptical  axes  at  a  given  x/L,  respectively 

Cy  Constant  used  In  k-e  turbulence  modeling,  generally  CH  •  0.09 

k  Turbulence  kinetic  energy  (k  -  (u^  +  v'1  +  w”‘i)/2) 

L  Total  body  length 

t  Mixing  length  parameter  (see  Equation  6) 

n,C  Coordinate  measured  normal  and  tangent  to  the  bodv  profile  In  the  yr 
plane,  respectively  (Fig,  I) 

r  Transverse  radius  of  curvature 

r0  Radial  distance  for  an  axlsymmetrlc  body 

rt  Transverse  radius  of  curvature  of  the  bodv 

llfl  Inviscld  flow  velocity  at  the  edge  of  the  boundary  layer 

I'l  Axial  component  of  the  Inviscld  flow  velocity  at  the  bodv  surface 

lt„  Free- st  ream  velocity 

0X  Axial  component  of  the  inviscld  flow  velocity 
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11  x*  vn.  WG  Mean  velocity  components  In  the  x,  n,  and  0 
t~  —  /— 

/(uV.  Av')\  A*')'  Turbulent  fluctuations  in  the  x. 


directions,  respectively 


n,  and  B  direction,  respectively 


Reynolds  stresses 


Turbulence  structurnl  parameter 
x,y,z  Orthogonal  body  coordinates  (Fig.  1) 

Yt  Interraittancy  Factor,  Tt  ”  (1  +  5.5  (— )b]*' 

fi 

6a,  6b  Boundary-layer  thicknesses  along  the  a  and  b  axes,  respectively 
ir  Boundary-layer  thickness  measured  In  the  n  direction 

ip*  Planar  displacement  thickness  (see  Equation  5) 

{995  Boundary- layer  thickness  at  which  ux/Ux  -  0.995 
c  Turbulent  dissipation  (see  Equation  3) 

0  Angular  coordinate  measure  In  the  y-z  plane 


0  Angle  between  the  x  and  0  coord lntes 

A*-Xu*/v  Generalized  Van  Driest  damping  factors 

X  Length  parameter  related  to  Van  Driest  damping  factors 

v  Kinematic  viscosity  of  the  fluid 

vt/(/k t)  Non-dimensional  turbulence  structure  parameter 

vc  Eddy  viscosity  (see  Equations  1  and  6) 

P  Mass  density  of  the  fluid 

xu  Wall  shear  stress 

1.  INTRODUCTION 

The  turbulent  flow  characteristics  in  a  ship's  stern  region  are  of  particular  Interest  to 
naval  architects  since  many  single-screw  propellers  operate  Inside  thick  boundary  layers. 
An  accurate  prediction  of  velocity  inflow  to  the  propeller  is  essential  to  meet  the  ever- 
increasing  demand  for  Improved  wake-adapted  propeller  performance.  The  integral  and/or 
differential  methods  used  to  predict  turbulent  boundary- layer  characteristics  agree  well 
with  experimental  data  in  thin  boundary  layers.  However,  predictions  from  these  methods 
deviate  from  measured  values  as  the  boundary  layer  thickens.  Often,  the  boundary  layer 
Is  thicker  than  the  transverse  dimension  of  the  model  in  the  stern  region,  nullifying  the 
haslc  assumption  that  the  boundary  layer  is  small  when  compared  with  model  dimensions.  It 
is  hoped  that  some  insight  Into  the  flow  in  thlB  complex  region  can  be  gained  through 
detailed  measurements  of  the  turbulent  houndary- layer  characteristics  in  the  stern  region 
of  simple  geometries.  To  minimize  the  variables  associated  with  three  dimensionality, 
models  with  2:1  and  3:1  elliptic  cross  sections  were  designed  to  have  the  same  longitudi¬ 
nal  distributions  of  cross-sectional  area  ns  one  of  the  axlsymmetrlc  models  tested 
[Axisymmetrlc  Model  11  in  Reference  1.  I'revlons  Investigations  at  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center  have  produced  comprehensive  sets  of  experimen¬ 
tal  data  for  pressure,  shear  stress,  velocity,  and  turbulence  across  the  sterns  of  three 
axlsymmetrlc  models  [l,  2  J ,  and  two  three-dimensional  models  having  3:1  and  2:1  elliptic 
transverse  cross  sections  [3,  4).  These  experimental  data  reveal  three  distinctive 
features  of  the  thick  stern  boundary- layer  flows:  (1)  strong  vlscous/invlscld  Interact¬ 
ion,  (2)  significant  cross-stream  pressure  variation,  and  (3)  special  mixing  length/eddy 
viscosity  with  a  distinct  "zonal"  turbulence  structure.  The  data  for  simple  three- 


u  v  ,  u  w 

-iT77(2k) 
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dimensional  models  having  1:1,  2:1,  and  3:1  elliptic  transverse  cross  sections  [1,3,4] 
are  used  to  show  the  special  turbulence  characteristics  of  the  stern  flows* 

In  the  following  sections,  the  models  and  experimental  techniques  are  given  in 
detail.  The  measurements  of  mean  velocities,  turbulence  Intensities,  and  Reynolds 
stresses  are  analyzed  to  obtain  eddy  viscosity  and  mixing  length.  The  measured  turbulent 
kinetic  energy  k,  and  the  derived  eddy  viscosity  vt,  the  mixing  length  t,  and  the  tur* 
bulent  dissipation  c  are  presented  to  demonstrate  that  the  thick  turbulent  stern  flows 
possess  a  distinct  "zonal"  turbulence  structure. 

2.  MODELS  AND  EXPERIMENTAL  TECHNIQUES 

The  detailed  offsets  of  Axisymmetric  Model  1  are  given  in  Table  1.  The  bow-entrance 
length/diameter  ratio,  parallel  middle  body  length/diameter  and  afterbody  length/diameter 
ratio  of  this  model  are  1.82,  4.85,  and  4.31,  respectively.  The  total  length  of  the  model 
is  3.066-m  (10.06-ft.).  The  diameter  of  the  parallel  middle  body  (rRax)  is  13.97-cm 
(0.458-ft).  The  longitudinal  distributions  of  the  transverse  cross  sectional  areas  of  the 
2:1  and  3:1  elliptic  models  are  equal  to  those  of  the  Axisymmetric  Model  1.  Thus,  the 
three  models  have  the  same  volume  and  longitudinal  distribution  of  buoyancy.  However,  the 
2:1  and  3:1  elliptic  models  have  9!  and  23X  more  wetted  surface  area,  respectively,  than 
the  axisymmetric  model.  A  schematic  of  the  three  sterns,  together  with  the  flow  measure¬ 
ment  stations,  is  shown  in  Fig.  1. 

The  experimental  investigation  was  conducted  in  the  DTNSRDC  Anecholc  Flow  Facility. 
The  wind  tunnel  has  a  closed  Jet  test  section  that  is  2.4-m  (8-ft)  square  and  4.19-m 
(13.75-ft)  long.  The  corners  have  fillets  which  are  carried  through  the  contraction.  The 
test  section  is  followed  by  an  acoustically  lined  large  chamber,  7.17-m  (23.5-ft)  long. 

It  was  found  previously  by  Huang  et  al.  [1]  that  the  ambient  free-streaa  turbulence 
levels  are  0.075,  0.090,  and  0.100,  and  0.12  ~  0.15  for  free-streaa  velocities  U0  of 
24.4,  30.5,  38.1,  and  45.7  m/s,  respectively.  Integration  of  the  measured  noise  spectrum 
levels  in  the  test  section  from  10  to  10,000  Hz  indicated  that  the  typical  background 
acoustic  noise  levels  at  30.5  m/a  were  about  93  dB  re  0.0002  dyne/cm^  (0.0002  Pa). 

These  levels  of  ambient  turbulence  and  acoustic  noise  were  considered  low  enough  so  as 
not  to  unfavorably  affect  the  measurement  of  boundary- layer  characteristics.  The 
maximum  air  speed  that  can  be  achieved  is  61  m/s  (200  ft/sec);  in  the  present 
experiments,  the  wind  tunnel  velocity  was  held  constant  30.48  m/s. 

The  models  were  supported  either  by  three  streamlined  guide  wires  with  one 
streamlined  strut  located  downstream  of  the  model  or  by  two  streamlined  struts  separated 
by  one-third  of  the  model  length.  The  disturbances  generated  by  the  supporting 
struts  were  within  the  region  below  the  horizontal  centerplane.  All  the  measurements 
were  made  in  each  model's  vertical  centerplane  along  the  upper  meridian,  where  there 
was  little  extraneous  effect  from  the  supporting  struts.  The  aft  one-half  of  the  model 
length  protruded  beyond  the  closed-let  section  into  the  open-jet  section.  The  ambient 
static  pressure  coefficients  across  and  along  the  entire  open-jet  chamber  (7.2  m  x  7.2  m 
x  6.4  m)  were  found  to  vary  by  less  than  0.3X  of  the  dynamic  pressure.  Tunnel  blockage 
and  longitudinal  pressure  gradient  effects  along  the  tunnel  length  were  almost  completely 
removed  by  testing  the  afterbody  in  the  open-jet  section.  The  two  three-dimensional 
models  are  shown  in  the  Anechoic  Wind  Tunnel  Facility  in  Figure  2.  The  location  of 
the  boundary- layer  transition  from  laminar  to  turbulent  flow  was  artificially  Induced  by 
a  0.61-mm  (0 . 024- in )-dlam  trip  wire  located  at  x/L  *  0.05.  Flow  visualization  techniques 
using  oil  dots  identified  a  small  region  of  flow  separation  on  the  two  three-dimensional 
sterns.  The  turbulent  flow  separation  on  the  3:1  model  was  found  to  begin  at  x/L  -  0.90 
and  to  reattach  at  x/L  »  0.98.  The  maximum  lateral  extent  of  the  separation  region 
occurred  at  x/L.  -  0.94  and  covered  a  region  of  80“<9  <90“.  The  measured  boundary  layer 
cross  flow  was  found  to  be  quite  large  (w/UQ  3  0.1)  at  9  *  80°  upstream  of  the 
separation.  The  smaller  separation  region  on  the  2:1  model  began  at  x/L  -  0.94  and 
reattached  at  x/L  »  0.96.  The  lateral  extent  of  the  separation  region  at  x/L  - 
0.95  was  around  83“<6<90“.  No  separation  was  found  on  the  stern  of  Axisymmetric 
Model  1 • 

The  mean  axial  and  radial  velocities  and  the  turbulence  intensities  for  the 
Fevnolds  stress  calculations  were  measured  by  a  TSI,  Inc.,  Model  1241-20  "X"  type 
hot-film  probe.  The  probe  elements  are  0.05  mm  (0.002  In.)  in  diameter  with  a 
sensing  length  of  1.0  rmn  (0.04  in.).  The  spacing  between  the  two  cross  elements 
is  1.0  mm  (0.04  in.).  A  two-channel  hot-wire  and  hot-film  anemometer  with  linearizers 
was  used  to  monitor  the  response  of  the  hot-film  probe.  A  temperature- compensat Ing 
sensor  ( probe )  was  used  with  each  hot-film  element  to  regulate  the  operating 
temperature  of  the  sensor  with  changes  in  air  temperature.  The  "X"  hot-film  and 
its  temperature-compensated  sensor  were  calibrated  together  through  the  expected 


air  temperature  range  and  so, -oiled  with  their  Individual  linearization  polynomial 
coefficients  at  the  factory# 

For  reliable  measurements,  the  frequency  response  of  the  anemometer  system  is 
claimed  by  the  manufacturer  to  be  0-100  kHz.  Calibration  of  the  "X"  hot-film  probe  was 
made  immediately  before  and  after  each  set  of  measurements.  The  maximum  deviation  of  the 
hot-film  output  from  the  mean  linear  response  curve  was  held  to  within  0.5Z  of  the  free- 
stream  velocity.  In  addition,  the  standard  deviation  of  the  hot-film  measurements  of  the 
free-stream  velocity  at  45.72  m/s  (150  ft/sec)  at  the  reference  location  were 
recorded  before  and  after  an  experimental  series  and  were  held  to  within  0.  5X.  An  esti¬ 
mate  was  made  of  the  cross-flow  velocity  by  yawing  the  "X"  hot-film  probe  at  the 
reference  free-ntream  location.  It  was  found  that  the  measured  cross-flow  velocities 
were  about  1Z  of  the  free-stream  velocity.  At  10  randomly  selected  points  across  the 
thick  boundary  layer,  the  streamwise  and  cross-flow  velocity  components  were  measured 
at  five  different  times.  The  standard  deviations  of  the  streamwise  and  cross-flow  velo¬ 
city  components  were  found  to  be  less  than  1Z  and  2Z,  respectively. 

The  linearized  signals  were  fed  into  a  Time/Data  Model  1923-C  real-time  analyzer. 
Both  channels  of  the  analog  signal  were  digitized  at  a  rate  of  128  points  per  second  for 
8  seconds.  These  data  were  immediately  analyzed  by  a  computer  to  obtain  the  individual 
components  of  mean  velocity,  turbulence  fluctuation,  and  Reynolds  stress  on  a  real-time 
basis. 


A  traversing  system  with  a  streamlined  strut  was  mounted  on  a  guide  plate  that  per¬ 
mitted  the  traverse  to  be  locked  into  various  stationary  positions  parallel  to  the  longi¬ 
tudinal  model  axis. 

All  data  are  presented  in  the  coordinate  system  used  to  experimentally  measure  the 
boundary- layer  flows.  The  coordinate  system,  denoted  by  x-n-0,  is  given  in  Fig.  1.  The 
axial  coordinate  x  Is  measured  from  the  nose  of  the  body  and  passes  through  the  center  of 
the  elliptic  profile.  The  coordinates  n  and  8  are  defined  along  an  axial  cut  normal  to 
the  x  axis,  l.e.,  in  the  yz  plane.  The  normal  component  n  is  measured  from  the  model  sur¬ 
face  and  Is  normal  to  the  elliptic  surface  but  not  to  the  model  as  a  whole.  The  angular 
coordinate  0  Is  defined  as  the  angle,  In  degrees,  measured  from  the  z-minor  axis  to  the 
line  joining  the  surface  offset  and  elliptic  center. 

Mean  velocity  and  turbulence  measurements  were  taken  with  an  ”X“  hot-film  sensor 
which  was  stepped  away  from  the  body  in  the  n  direction.  Measurements  of  velocity  in  the 
axial  x  and  normal  n  directions,  ux  and  vn,  respectively,  were  taken  with  the  probe  ele¬ 
ments  aligned  vertically.  The  sensor  elements  were  rotated  90*  to  the  horizontal  posi¬ 
tion  to  measure  the  mean  velocity  wg  in  the  8  direction.  An  on-line  computer  was  used  to 
collect  data  at  a  sample  rate  of  1024  data  Values  in  8  sec.  The  root-mean  square  values 
of  turbulence  velocity  were  recorded  at  each  probe  position,  and  the  eddy  viscosity  and 
mixing  length  values  were  computed  from  the  measured  Reynolds  stresses  and  the  measured 
mean  velocity  profiles. 


The  distribution  of  the  Reynolds  stresses  -u^v”,  -u^w’’,  (u^)  ,  (v')Z,  and  (w”  )* 
represent  the  turbulence  characteristics  In  the  thick  boundary  layer.  The  mean-square 

turbulent  velocity  fluctuations  (u* )'  In  the  axial  direction  and  (v*  In  the  n  direc¬ 
tion,  and  the  Reynolds  stress  -u  v  were  measured  with  the  "X”  hot-film  probe  elements 
aligned  vertically.  The  probe  elements  were  rotated  90"  to  the  horizontal  position  to 


measure  both  the  turhulent 


f luctuat Ion 


in  the  6  direction  and  the  Reynolds  stress 


-  u 


w  , 


I. Inear  Interpolation  was  used  to  approximate  (w")?  and  -u"w"  at  the  same  off- body  posi¬ 
tions  as  the  data  measured  In  the  vertical  direction. 

1.  MKAStlKKI)  AND  DKRIVKI)  TIIKBIH.KNT  CHARACTERISTICS  OF  THICK  STERN  FLOWS 

The  measured  and  derived  turhulent  characteristics  relevant  to  the  analysts  of  thick 
stern  boundary-layer  flows  are  presented  below. 


3.1  Measured  Reynolds  Stresses 


The  typical  turbulence  characteristics  in  the  thick  boundary  layer  can  be 

represented  by  the  distributions  of  Reynolds  stresses,  namely,  -u'* v"”  ,  where 

u  ,  v  ,  and  w  are  the  turbulent  velocity  fluctuations  in  the  axial,  radial,  and  azimuthal 
directions,  respectively.  Figure  3  shows  the  measured  distribution  of  Reynolds  stress 

-u  v  /U0^  and  three  components  of  turbulence  intensity  at  several  axial  locations  along 
the  1:1  afterbody.  In  general,  for  a  given  location,  the  intensity  of  the  axial  tur¬ 
bulent”  velocity  component  has  the  highest  value  and  the  intensity  of  the  radial  com¬ 
ponent  has  the  smallest  value.  The  degree  of  anisotropy  decreases  as  the  stern 
boundary  layer  becomes  thicker.  Furthermore,  the  increased  boundary- layer  thickness  is 
accompanied  by  a  reduction  of  turbulence  intensities  in  the  inner  region.  The 
variation  along  the  body  of  the  radial  location  of  the  maximum  values  of  the  trea¬ 
sured  Reynolds  stress  -u*v*  /U0^  layer  is  small.  The  spatial  resolution  of  the 
cross-wire  probe  is  not  fine  enough  to  measure  the  Reynolds  stress  distribution  in  the 
inner  region  when  the  boundary  layer  is  thin.  As  the  stern  boundary  layer  increases 
in  thickness,  the  location  of  maximum  Reynolds  stress  moves  away  from  the  wall  (Figure 

3).  The  values  of  Reynolds  stress  -u  v  decrease  quickly  from  the  maximum  value  to  zero 
at  the  edge  of  the  boundary  layer.  As  shown  in  Figure  3,  the  shape  of  the  Reynolds 
stress  distribution  curves  in  the  outer  region  is  quite  similar  for  all  the  thick  boun¬ 
dary  layers.  It  is  interesting  to  note  that  the  shapes  of  the  measured  Reynolds  stresses 
In  the  inner  regions  of  the  near  wake  at  x/L  *  1.057  and  1.182  experience  a  drastic 
reduction  in  magnitude. 

3.2  Measured  Turbulence  Kinetic  Energy  k 

The  measured  turbulence  kinetic  energy  k  *  ( u  ^  +  v'1  ♦  w**)/2  is  shown  in  Figure  4.  The 

thin  boundary  layer  can  be  seen  in  the  region  of  x/L  <  0.8  of  the  1:1  model  and  as  the 
large  flat  portions  of  the  2:1  and  3:1  models,  0  <  60°.  The  maximum  value  of  k 
measured  In  this  region  occurs  very  close  to  the  wall  (less  than  10Z  of  the  boundary- 
layer  thickness)  and  the  measured  value  of  k  decreases  almost  linearly  towards  the  edge 
of  the  boundary  layer.  However,  the  measured  k  in  the  Inner  region  of  the  thick 
boundary  layer  has  a  constant  value  for  almost  60Z  of  the  boundary- layer  thickness.  The 
distinct  characteristics  of  the  measured  values  of  k  in  the  thick  boundary  layer  are 
found  In  Fig.  5  for  x/L  >0.9  for  the  1:1  model  and  In  the  corner  region  of  the  major 
axis  (75°<0  <90°)  of  the  2:1  and  3:1  models. 

3.3  Derived  Eddv  Viscosity,  Mixing  Length,  Turbulence  Kinetic  Energy  Dissipation 

The  values  of  eddy  viscosity  and  mixing  length  are  not  measured  directly  but  are  obtained, 

as  In  previous  studies  [ 1-4]  from  the  measured  values  of  the  Reynolds  stress  -u  v  and  the 
mean  velocity  gradient  9ux/9n.  The  definitions  used  to  compute  these  quantities  are 


When  the  values  of  wq/Ux  are  less  than  0.1  and  0,  the  angle  between  the  coordinates 
x  and  O  (Fig.  1)  Is  90°  for  the  present  measurements.  Equation  (1)  may  be  approximated 
by 


(2) 


A  spline  curve  is  used  to  fair  the  experimental  data  before  the  velocity  gradient  is 
obtained  numerical  1 v.  The  dissipation  of  the  turbulence  kinetic  energy  is 
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It  ts  necessary  to  measure  each  of  the  above  nine  mean-square  derivatives  to  determine 
accurately  the  value  of  e.  Because  of  the  practical  limitations  of  hot-wire  techniques  and 
the  tremendous  time  required  to  measure  even  a  few  of  the  nine  terms,  a  direct  measurement 
of  the  value  of  e  was  not  attempted.  However,  bv  dimensional  analysis  of  k  and  e,  one 
finds 


where  the  constant  *  0,09  is  recommended  hv  Launder  et  al .  [5, 6),  and  is  generally 

accepted  in  the  community  of  k-e  turbulence  modeling.  Since  the  values  of  k,-u  v  ,  and 
ux  (nence  3ux/3n)  were  determined  hv  direct  measurements  and  the  values  of  vt  were  deri¬ 
ved  from  Equation  (2),  the  values  of  e  can  then  be  derived  from  Equation  (A). 

4.  COMPARISON  OP  SIMPLE  THIN  BOUNDARY  LAYER  TURBULENCE  MODELS  WITH  THE  MEASURED  DATA 

The  simple  mixing  length  model  of  Bradshaw  et  al.  f  7 ]  and  the  simple  eddv  viscosity 
model  of  Cebeci  and  Smith  [8]  are  compared  with  the  present  measured  data  from  the  thick 
stern  boundary  layer. 

4.1.  Comparison  of  the  Bradshaw  et  al.  [7)  Mixing  Length  Model  With  Present  Data 

When  the  past  history  of  boundary  laver  characteristics  is  important,  Bradshaw  et  al. 
(7)  argue  that  the  turbulence  energy  equation  can  be  used  to  model  the  "memory”  effect.  In 
order  to  determine  the  rate  of  change  of  the  shear  stress  along  a  mean  streamline,  three 

assumptions  are  made:  namelv,  (1)  the  local  Reynolds  stress  -u~ v' /(2k)  «  0.15;  (2)  the 
dissipation  rate  is  determined  by  the  local  Reynolds  stress  and  a  length  scale  depending 
on  n/6r,  and;  (3)  the  energy  diffusion  is  directly  proportional  to  the  local  Reynolds 
stress  with  a  factor  depending  on  the  maximum  value  of  the  Reynolds  stress.  On  the  basis 
of  thin  boundary- laver  data,  two  empirical  functions  for  the  last  two  assumptions  were 
proposed  by  Bradshaw  et  al.  (7),  The  first  assumption  t/6  «  (nj6r)  can  be  evaluated  by 
the  present  thick  boundary  laver  data  and  is  shown  In  Fig.  5.  This  assumption  is  found 
not  to  be  applicable  to  thick  stern  flows.  The  measured  values  of  £/  6  r  are  much  smaller 
(sometimes  hv  a  factor  of  3)  than  those  of  the  thin  boundary  laver. 

It  is  important  to  note  that  the  boundar v 1  aver  thickness  of  a  typical  body  increases 
drastically  at  the  stern.  Most  of  the  rapid  change  takes  place  within  a  streamwlse 
distance  of  a  few  houndarv- laver  thicknesses.  Most  of  the  empirical  functions  for 
solving  the  turbulence  enerev  equation  will  undergo  rapid  changes.  The  one  known  func¬ 
tion  Is  the  empirical  function  for  mixing  length,  as  shown  in  Fig.  5,  for  defining  the 
value  of  £/<Sr  as  a  function  of  n/<5r.  Therefore,  the  empirical  functions  used  by  Bradshaw 
et  al.  [7 |  need  to  be  revised  for  the  thick  stern  boundary  lavers. 


4.2  Comparison  of  the  Cebeci  and  Smith  (8)  Eddv  Viscosity  Model  With  Present  Data 


In  most  simple  turbulence  models,  the  basic  assumption  made  In  the  differential 
methods  for  calculating  turbulent  boundary  layers  is  that  the  Reynolds  stress  is  modeled  hv 
the  product  of  the  mean  velocity  gradient  and  the  eddv  viscosity  or  the  square  of  the  mean 
velocity  gradient  and  th-  mixing  length.  Both  the  eddv  viscosity  and  the  mixing  length 
are  uniquely  related  to  the  boundary- 1  aver  thickness  parameter  at  the  local  station.  As 
long  as  the  boundary  laver  is  thin  and  the  change  in  boundary- laver  properties  due  to  the 
pressure  gradient,  and  bodv  geometry  is  gradual,  this  simple  turbulence  model  is  known  to 
he  sat isf actorv ,  (e.p.  fleheci  and  Smith  |H|).  The  exper Lmenta 1 1 v-determined  distributions 
of  eddv  viscosity,  are  shown  In  Fig,  b  for  the  three  sterns,  where  6p*  is 

the  displacement  thickness,  defined  as 
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v„  -  0.0168  US  «p*  y tr 

A  ■  26v  (Tw/p)_1*  -  26v/u*,  Van  Driest  damping  factor 

Ytr  "  [1  +  5.5  (n/df)6]'1 

6r  -  6995,  boundary  layer  thickness  at  which  ux/Ux  ■  0.995 

tw  -  wall  shear  stress 

As  shown  In  Figure  6a,  the  measured  distribut ions  of  eddv  viscosity  agree  reasonably 
well  with  the  eddy  viscosity  of  Cebeci  and  Smith  [8]  when  the  boundary  layers  are  thin. 
However,  as  the  stern  boundary  layer  thickens,  the  measured  values  of  vt/U5$p*  in  the  thick 
stern  boundary  layer  shown  in  Fig.  6  are  significantly  (as  much  as  one  order  of  magnitude) 
smaller  than  the  values  for  thin  boundary  layers  given  by  the  Cebeci  and  Smith  model. 
Therefore,  this  commonly-used,  simple  turbulence  model  of  Cebeci  and  Smith  for  thin 
boundary  layers  is  not  suitable  for  thick  boundary  layers. 

5.  -ZONAL**  TURBULENT  SIMILARITY  PROPERTIES  OF  THICK  STERN  FLOWS 

The  solutions  of  the  t ime* dependent ,  three-dimensional,  Navier-Stokes  equations  with 
adequate  resolution  of  time  and  spatial  scales  by  today's  computers,  are  limited  to  lamiftar 
flows  at  very  low  Reynolds  numbers.  The  time-averaged  Navier-Stokes  equations  are  often 
used  to  compute  the  mean  flow  quantities  at  high  Reynolds  numbers,  but  the  time-dependent 
solutions  are  then  sacrificed.  The  simplified  time- averaged  Reynolds  equations  require  a 
set  of  closure  relationships  for  Reynolds  stresses.  The  closure  relationship  reduces  the 
number  of  unknowns  to  the  number  of  equations.  The  common  approach  to  closure  is  to 
define  an  eddy  viscosity  (an  artifice)  for  turbulent  flows  in  the  same  form  as  the  laminar 
viscosity  in  the  stress  tensor  for  the  Newtonian  flows.  The  specification  of  the  eddy 
viscosity  in  terms  of  algebraic  or  differential  equations  has  been  the  major  area  of 
turbulence  research  for  the  past  20  years.  Some  of  the  turbulence  models  are  sue* 
cessful  for  thin  boundary  layers  and  for  simple  flows,  but  none  of  the  closure  rela¬ 
tionships  is  universally  valid  for  all  complex  flows.  Therefore,  Kline  (10)  in  his 
summary  of  the  1980-81  Stanford  Conference  emphasized  the  importance  of  developing  "zonal 
models"  for  various  structural  flow  zones.  Thick  turbulent  stern  flows  possess  a  special 
similarity  flow  structure  which  differs  from  that  of  thin  boundary  layers.  These 
distinct  "zonal"  characteristics  of  turbulent  stern  flows  are  discussed  below. 

5,1.  Similarity  Property  of  Mixing  I-ength  for  Turbulent  Stern  Flows 

For  three  axisymmetric  and  two  simple  three-d intensions  1  turbulent  stern  flows,  Huang 
et  al.  (1,2,3,41  obtained  a  similarity  property  of  the  measured  mixing  length.  The  mixing 
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length  was  found  to  he  proportional  to  the  square  root  of  the  entire  turbulence  area 
between  the  body  surface  and  the  edge  of  the  boundary  layer.  Although  the  apparent  boun¬ 
dary  layer  thickens  rapldlv  toward  the  stern,  this  turbulent  area  undergoes  a  relatively 
small  change  In  the  stern  region.  Therefore,  the  mixing  length  does  oof  change  appre¬ 
ciably  in  the  stern  region.  As  shown  in  Figures  3  and  4,  the  measured  values  of  Reynolds 
stresses  and  turbulent  kinetic  energy  k  in  the  thick  stern  boundary  layer  decrease  from  a 
relatively  large  value  at  60  percent  of  the  boundary  layer  thickness  to  zero  at  the  outer 
edge  of  the  boundary  layer.  Contrary  to  the  thin  boundary  layer  data,  the  measured 
values  of  k  in  the  thick  boundary  laver  remain  almost  constant  within  60  percent  of  the 
boundary- layer  thickness  (Fig.  4),  An  approximate  correlation  of  the  measured  contours 

of  non-dimensional  turbulence  kinetic  energy  /k/0o  and  the  measured  contours  of  the 
boundary- laver  thickness  for  the  2:1  model  is  shown  in  Fig.  7.  The  measured  contour  of 

0.66r  i s  very  close  to  the  measured  contour  of  /k/0o  =  0,04,  Within  the  contour  of  0*66r 
and  the  body  surface,  the  measured  values  of  /k/U0  are  found  to  vary  from  0.035  to  0.055 
for  all  the  data  taken  for  the  three  models.  The  measured  values  of  /k/UG  reduce  rapidly 

to  about  0.01  at  the  edge  of  the  boundary  laver. 

For  axisymmetric  and  simple  three-dimensional  elliptic  sterns,  the  gross  turbulence 
area  A(x)/w  =  [  (a+0. 66a) (b+0. 66 ^,)  -  (a+€a) (b+e^,)  ] ,  may  be  more  relevant  to  the  mixing 
length  parameter,  where  ea  and  are  the  effective  thickness  of  the  separation  bubble  (to 
allow  for  low  turbulence  mixing)  in  the  direction  of  the  major  and  minor  axes,  a  and  b, 
respectively,  of  the  elliptical  cross-section,  and  6a  and  <5^  are  the  boundary-layer 
thicknesses  along  the  a  and  b  axes.  A  small  separation  bubble  (covering  0. 90<  x/L<0. 98  for 
the  3:1  model  and  0. 92< x/L<0. 96  for  the  2:1  model)  was  detected  in  the  corner  region  of 
the  major  axis,  and  the  values  of  ea  are  found  to  be  negligibly  small.  The  flow  does  not 
separate  from  the  entire  flat  portion  of  the  two  elliptic  sterns  along  the  minor  axis  (e^  ■ 
0).  For  the  axisymmetric  stern  a*b*r0,  ^a^b^r*  and  ea  *  *  0  (the  flow  remained 

attached  over  the  entire  body  for  axisymmetric  models  tested).  The  similarity  property  of 

the  mixing  length  i  for  the  turbulent  stern  flows  can  be  evaluated  by  plotting  i//A/w 
versus  n/6r,  where  n  Is  the  distance  normal  to  the  body  surface  and  the  longitudinal  axis 
6r  is  the  boundary-layer  thickness  in  the  n-direction,  taken  as  the  value  of  n  at  which 
ux/Ux  -  0.995,  ux  is  the  viscous  axial  velocity  and  Ux  is  the  inviscld  axial  velocity. 

At  the  location  of  6r,  the  measured  value  of  /k/UQ  Is  about  0.01. 

The  normalized  mixing  length  distributions  for  the  three  axisymmetric  bodies 
investigated  by  Huang  et  al.  (1,2)  are  shown  in  Figures  8a  through  8c;  Figures  8d  through 
8g  show  the  distribution  for  the  2:1  model  (4);  and  Figures  8h  through  8j  show  the 
distributions  for  the  3:1  model  (3).  As  shown  in  Figure  8,  the  measured  values  of  the 
normalized  mixing  length  are  fairly  similar  over  the  entire  thick  boundary  layer  region  of 
all  the  sterns.  This  similarity  property  remains  valid  even  for  the  mild  separation 
region  of  the  elliptic  models.  The  curve  shown  in  Figure  8  is  the  average  line  of  the  data 
and  can  be  represented  by 


0. 1 69  — ! —  exp 
*r 


-  —  .21  /n\ 3 

*>  5r  in  '  5rJ 


(O 


It  Is  assumed  that  Equation  (7)  is  approximately  valid  for  an  arbitrary  three- 
dimensional  stern.  The  boundary-layer  thickness  <5r  (and  hence  0.66r)  can  be  obtained  from 
most  boundary- layer  computation  procedures.  Therefore,  a  simple  numerical  evaluation 
can  be  performed  to  compute  the  effective  turbulence  area  A /it  for  general  three- 
dimensional  stern  flows. 

Tt  is  important  to  note  that  when  the  boundary- layer  thickness  of  the  axisymmetric 
stern  equals  23  percent  of  the  transverse  radius  of  curvature,  rt=a“b=r0=6 r/0.23 ,  Equation 
(7)  becomes 

t„  -  0.40  n  exp  |  -  ii 

*>  5r  10 

This  relationship  very  closely  resembles  the  empirical  mixing  length  curve  proposed  by 
Bradshaw  et  al .  [7|.  As  shown  In  Fig.  5a,  Equation  (8)  obtains  a  maximum  £0/6r  of 

0,096  at  nMr  =  0.5.  Since  Equation  (8)  and  the  Bradshaw  et  al.  empirical  mixing 


8 


length  curve  are  In  good  agreement,  Equation  (8)  will  be  assumed  valid  for  all  thin  boun¬ 
dary  layers  with  the  boundary  layer  thicknesses  less  than  21  percent  of  the  transverse 
radii  of  curvature  of  the  bodies. 


Although  Equations  (7)  and  (8)  specify  the  mixing  length  distribution  for  the  entire 
thick  and  thin  boundary  layers,  respectively,  the  classical  mixing  length  (consistent  with 
the  logarithmic  velocity  laws)  Is  more  appropriate  for  the  inner  wall  region.  For 
example,  the  generalized  mixing  length  proposed  by  Granville  111), 

*1  *  j=±.  I  1  *  exp  (  -  y*/X*)  ]  (9) 


X*  -  2b  +  ib  X  1113  I  +  2.71  X  in'*;  i  \ 


can  be  used,  where  rt  Is  the  transverse  radius  of  curvature,  y  is  the  distance  normal  to 

the  wall,  X*  is  the  generalized  Van  Driest  damping  factor,  X*“Xu*/v,  >'s*“bu»/v,  u*  -  / T„/p 
Tw  is  the  wall  shear  stress,  v  is  the  kinematic  viscosity,  p  is  the  mass  density 
of  the  fluid,  and  X  is  the  length  factor.  Equation  (9)  reduces  to  the  classical 
mixing  length  for  the  thin  two-dimensional  boundary  layer  when  rt  ♦  “,  e.g.. 


0. Ay  [1  -  exp(-y*/26)) 


The  value  of  Increases  almost  linearly  with  increasing  y  for  large  values  of  y.  In  the 
outer  boundary  layer,  the  maximum  value  of  the  mixing  length  is  limited  by  Equation  (7) 
for  the  thick  boundary  layer  and  by  Equation  (8)  for  the  thin  boundary  layer.  At  the 
intersection  of  the  inner  and  the  outer  regions  (n-nc  shown  in  Fig.  9),  the  condition  of 
l0-tj  must  be  applied.  The  shaded  line  shown  in  Fig.  9  is  the  proposed  mixing  length 
distribution  across  both  the  thick  and  thin  boundary  layers,  and  may  be  summarized  as 
follows, 

{Equation  (7)  for  thick  boundary  layer,  Sr>0.23  rt  (10) 

for  outer  region. 

Equation  (8)  for  thin  boundary  layer,  *rc0.23  rt  n  >  nc 


tj  “  Equation  (9) 


for  inner  region, 
n  <  n. 


and  £0  -  at  n  “  nc 


The  boundary  layer  is  classified  as  thick  when  the  boundary- layer  thickness  ts  equal  to  or 
greater  than  23  percent  of  the  transverse  radius  of  the  curvature  of  the  body.  For  a  small 
separation  bubble  at  the  stern.  Equations  (7)  and  (8)  can  still  be  applied  outside  of  the 
separation  bubble.  However,  Equation  (9)  gives  tjO  at  separation,  due  to  ue  “  0, 
and,  therefore,  becomes  an  unnecessary  inner  region. 

Existing  thin  boundary- layer  methods  are  well  developed  and  can  be  applied  from  the 
forward  portion  of  the  body  to  the  location  where  the  boundary- layer  thickness  Increases  to 
23  percent  of  the  radius  of  the  transverse  curvature  of  the  body.  Further  downstream,  the 
proposed  similarity  property  of  the  thick  stern  boundary  layer  given  by  Equation  (10)  can 
be  used  to  solve  the  turbulent  Btern  flows.  The  presently  proposed  direct  model  for 

calculating  £  (and  hence  vt  and  -u’ v* by  Equation  (10))  without  the  need  of  solving  a 
differential  equation  can  be  called  a  Zero-Equat  ion  Model. 

5.2  Measured  Turbulence  Structural  Parameters 

In  turbulence  modeling,  the  non-dimensional  turbulence  structural  parameters  -u^v”  / k 

and  vt/(/k i)  are  often  constant  In  the  turbulent  flow  field.  The  proportionality 
constants  obtained  for  the  near-wall  turbulence  are  assumed  to  be  universal  coefficients 


and  to  be  valid  for  general  turbulent  flows.  These  proportionality  constants  have 
been  derived  by  Launder,  et  al.  (5,6,12)  as  well  as  many  others.  A  brief  summary  of  the 
characteristics  of  the  near-wall  turbulence  will  be  repeated  and  the  proportionality 
constants  obtained  from  the  near-wall  turbulence  will  be  compared  with  the  measured  values 
for  the  entire  thick  stern  flows. 


Very  near  the  wall,  the  stress  li/u^  is  dominated  by  u  v  ;  and  the  turbulent  velocity 
fluctuations  (with  turbulence  fluctuations  u^In  the  primary  flow  direction  x  and  v  in  the 
primary  velocity-gradient  direction)  are  small,  while  the  turbulence  production  and  dissi¬ 
pation  are  large.  Therefore,  the  advection  term  and  the  diffusion- loss  term  in  the  tur¬ 
bulence  energy  equation  [7]  can  be  neglected,  e.g.. 


Iul!l  I 

3** 

reduces  Co  a  balance  of 


k  - 

the 


T  3ux  3  .  -  .  ..  .  lui'  3u/ 

-  — -  +  -  (  P  V  +  kv  ]  +  — l  — L 

p3n3n 

turbulence  energy  production  and  dissipation,  or 


1  3ux 
p  dn 


V  3u  1  3ul 
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Furthermore,  the  convective  acceleration  term  and  Che  longitudinal  pressure  gradient  are 
usually  much  smaller  than  the  normal  shear  stress  gradient.  Hence  the  mean  flow  momentum 
equation 


3 

3 


I  ux 


.  i  (  l£) 

p  v  ax' 


reduces  to 


-  0 


which  means  that 


_T_ 

P 


P 


or 


u* 


2 


(11) 


if  the  small  contribution  of  the  laminar  stress  is  neglected. 
Using  the  definition  of  mixing  length 


“  u* 


2 


the  approximated  velocity  gradient  near  the  wall  becomes 
3ux  m  u* 

3n  t 

The  approximated  turbulence  energy  dissipation  near  the  wall  can  be  obtained  as 

t  3uv  -r-T  3u- 

e  »  —  _ _  *  -  u  v  _ z  =  _ 

p  3n  3  n  t  (12) 


If  the  commonly  used  relationship  of  eddy  viscosity  (with  Cy  -  0,09)  is  applied,  a 

very  useful  result  Is  obtained,  e.g.. 


10 


-u  V 


Cu  k.2 

3  n 


Therefore,  one  of  the  proportionality  constants  of  the  turbulence  structural  parameter  can 
be  evaluated, 


fr  or 


0.15  (with  Cy  -  0.09) 


The  other  proportionality  constant  of  the  turbulence  structural  parameter  is 


/■u*2 


CM  k V2 


!/4  -  0.  548 


have  used  the  alternative  turbulence  structural  parameter 

[U*2  \  3/2 

k3/2  _  VcT/  .  -3/4 


The  typical  measured  distributions  of  -u  v  /(2k)  are  shown  in  Fig.  10.  The  present 

thick  stern  flow  data  shown  in  Fig.  10  indicate  that  the  measured  values  of  -u' v' /(2k)  are 
almost  constant  at  an  approximate  value  of  0.15  up  to  0.66r  and  decrease  towards  the  edge 
of  boundary  layer.  It  should  be  pointed  out  that  the  measured  values  of  k  Include  the 
contribution  of  free-stream  turbulence,  no  attempt  having  been  made  to  remove  the  free* 
stream  turbulence  from  the  measured  values  of  k.  The  measured  reduction  of  the  propor* 
tionality  constant  near  the  outer  edge  of  the  boundary  layer  is  in  part  caused  by  the 


larger  contribution  of  the  free-stream  turbulence  to  k  than  to  *u  v 


Nevertheless,  the 


measured  values  of  the  turbulence  structure  parameter  of  -u  v  /(2k)  are  quite  constant  and 

agree  well  with  the  proport  Iona  1 i tv  constant  of  -u  v  /(2k)  •  0.15  derived  from  the  near- 
wall  turbulence  (Equation  (13))  within  W>  percent  of  the  boundary  layer  thickness,  the 
region  where  the  effect  of  free-stream  turbulence  is  small. 

The  measured  distributions  of  vt/(/k l)  are  shown  in  Fig.  11.  The  proportionality 
constant  vt/(/k l) y  derived  from  the  near-wall  turbulence  (Equation  (14)),  is  0.548.  The 
present  thick  stern  flow  data  shown  in  Fig.  11  Indicate  that  the  measured  values  of  vt(/kl) 

are  nearly  constant  for  the  entire  boundary  laver.  The  measured  values  of  v^/(/ki)  for 
1:1  model  agree  very  well  with  the  proportionality  constant  of  0.548  derived  from  the 

near-wall  turbulence.  However,  the  measured  values  of  v^/(/kt)  for  the  2:1  and  3:1  models 
are  slightly  lower  than  0. 548  and  decrease  slightly  towards  the  edge  of  the  thick  boundary 

laver.  The  consistency  of  the  measured  values  of  vt/(/kt)  *  Cy ^  *  0.548  across  the 
entire  thick  boundary  laver  can  he  used  to  estimate  the  distribution  of  the  turbulent  kine¬ 
tic  energy,  k=(vt/t)^  /cu  once  the  distribution  of  the  values  of  1  and  hence  the  values 
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of  v£  are  known  across  the  thick  stern  flow.  Furthermore,  the  distribution  of  the  tur- 
bulence  energy  dissipation  can  be  estimated  by  £  *  Cy  =  vtVl^. 

6.  CONCLUSION 

Direct  measurements  of  Reynolds  stresses,  turbulence  kinetic  energy,  and  mean  velocity 
profiles  have  been  made  by  “X"  hot-film  anemometry  for  the  stern  region  of  simple  three- 
dimensional  models  having  1:1,  2:1,  and  3:1  elliptic  transverse  cross  sections.  The 
values  of  eddy  viscosity  and  mixing  length  were  obtained  from  the  measured  Reynolds 
stresses  and  the  mean  velocity  gradients.  The  following  conclusions  can  be  drawn  from  the 
detailed  analyses  of  the  directly  measured  experimental  data  on  the  turbulence  charac¬ 
teristics  of  thick  stern  flows. 


The  experimentally-determined  distributions  of  eddy  viscosity  and  mixing  length  in  the 
thick  stern  boundary  layers  were  found  to  be  smaller  than  values  which  have  been  proposed 
for  thin  boundary  layers.  A  special  similarity  property  of  mixing  length  for  turbulent 
stern  flows  has  been  found.  The  mixing  length  in  the  outer  region  of  the  thick  boundary 
layer(l0)  is  found  to  be  proportional  to  the  square-root  of  the  effective  turbulence  area 
A/s,  where  A  la  the  area  between  60  percent  of  the  boundary-layer  thickness  and  body's  sur¬ 
face.  The  mixing  length  distribution  i  in  the  outer  region  is  given  by 
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0.169  /  A/x 
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ir  >  0.23  rt 


r,  /  .on 

0.4  n  exp  I  -  —  —  - 
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n  _  J52  Jn  \ 

5r  30  Ur' 


for  thin  B.L. ,  i ^  0.23  rt 


where  n  is  the  distance  normal  to  the  surface  and  the  longitudinal  axis  and  (r  is  the 
boundary-layer  thickness  in  the  n-directlon,  taken  as  the  value  of  n  at  which  the  boundary 
layer  velocity  approaches  0.995  of  the  lnvlscid  velocity.  The  classical  mixing  length 
(t^),  consistent  with  the  logarithmic  velocity  laws  is  more  appropriate  for  the  inner 
region  (Equation  (9)),  where  the  conditions  i  -  ij  for  n<nc,  l  -  to  for  n>nc  and  t„  ■  ij 
at  n  *  nc  are  imposed,  for  the  application  of  the  revised  mixing  length  model. 


Across  60  percent  of  the  thick  stern  boundary  layer,  the  measured  turbulence  struc¬ 
tural  parameter,  -iT’v’ 7(2k)  was  found  to  be  approximately  0.15,  a  value  in  agreement  with 


-u’v*/(2k )  -  /c^/2  derived  from  the  near-wall  turbulence.  Another  measured  turbulence 
structural  parameter  vt/(/kt),  was  found  to  be  almost  constant  across  the  entire  thick 


stern  boundary  layer.  The  proportionality  constant  for  vt/(/ki),  given  by  C y  ^  *  0.548,  was 
derived  from  the  near-wall  turbulence  structure,  and  found  to  be  approximately  valid  for 
the  entire  stern  boundary  layer.  The  consistency  of  these  two  measured  turbulence 
structural  parameters  together  with  the  definition  of  vt.Cykz/e  can  ^  used  to  cstimate  the 


distribution  of  turbulence  energy,  k 
energy  dissipation,  e  -  vtVt^ 
l  and  hence  the  eddy  viscosity  (vt 
layer. 


and  the  distribution  of  turbulence 


C„^'*|^'*/i  as  long  as  the  distributions  of  mixing  length 
1  3ux/3n)  are  known  across  the  entire  thick  boundary 
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Figure  7  -  Representing  the  Mixing  Length  of  the  Stem  Flow  by  the  Square  Root 

of  Turbulence  Area 


Figure  9  -  Schematic  of  Mixing  Length  across  the  Boundary  Layer 


Figure  10  -  Measured  Distribution  of  Turbulent  Structure  Paraaeter 
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g— iwr.y.j 

0.  0370 

0.4076 

0.0456 

0.8765 

niKi;  iW^I 

0.0746 

0.0381 

0.4175 

0.8841 

0.0795 

0.  4274 

0.0955 

0.0645 

M— WTiWi<' jjJ 

0.4374 

0.0456 

Mil  i  1  1 

"Si 

0.  0407 

0.  4473 

0.0456 

0. 9144 

0.0211 

0.0944 

0.0414 

0.  4573 

0.0456 

0.9245 

0. 0109 

0.0421 

0.4672 

0.9344 

0.0166 

0.  1044 

0.  0427 

0.  4771 

0.0456 

0.9443 

0.0140 

0.  1093 

0.0432 

0.4871 

mmwim 

0.9613 

0.0122 

0.  1143 

0.0437 

0. 4970 

0.0456 

0.9563 

0.0108 

0.  1193 

0.0441 

0.5070 

0.9612 

0.0095 

0.1243 

0.0444 

0.5169 

0.0456 

0.9642 

0.0087 

0.  1292 

0.  0447 

0.5268 

0.0456 

0.9662 

0.  0081 

0.  1342 

0.0450 

0.  5368 

0.0456 

0.9682 

0.0076 

0. 1392 

0.0452 

0.5467 

jfsKjFTITi 

0.9692 

0.0074 

0.  1441 

0.0453 

0.5567 

0.0456 

0. 9702 

0.0072 

0. 1491 

0.0454 

0.5666 

0.0456 

0.9722 

0.0066 

0.  1541 

0.0455 

0.5765 

0.0456 

0.9732 

0.0066 

0. 1590 

0.0456 

0.5865 

0.0456 

0. 9751 

0.0083 

0. 1640 

0.0456 

0.5964 

0.0456 

0. 9771 

0.0062 

0. 1690 

0.0456 

0.0456 

0.9791 

0.0059 

0. 1740 

0.0456 

0.6188 

0.0450 

0.9811 

0. 0056 

0.1789 

0.0456 

0.6264 

0.0455 

0.9831 

0. 0053 

0.  1839 

0.6378 

0.0455 

0.9851 

SHiHnfMTiii 

0.  1889 

0.0456 

0.6454 

0.0455 

0.9871 

0.0048 

0.  1938 

0.0456 

0.6567 

0.0453 

0.9881 

0.0046 

0.  1988 

0.0456 

0.6681 

0.0452 

0.9901 

0.0043 

0.2087 

0.0456 

0.  6757 

0.0450 

jSHHtWMTTjj 

0.2187 

0.0456 

0. 6871 

0.0448 

0.9940 

0.0036 

0.2286 

0.0456 

0.6984 

0.0444 

0.0028 

0.2386 

0.0456 

0.  7060 

0.0441 

0.0019 

0.2485 

0.2584 

0.0456 

0.0456 

0.7174 
0.  7250 

0.0437 

0.0433 

0.0000 

TABI.F.  1  -  Offset  of  DTNSKDC  Ax  l.symni-tr  ic  Modrl  I 
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INITIAL  DISTRIBUTION 


ies 

1  WES 

1  U.S.  ARMY  TRANS  R&D 

Marine  Trans  Div 

3  ONR/432F  Whitehead, 

Lee,  Reischman 

I  ONR/Bosuon 

1  ONR/Chicago 

1  ONR/New  York 

I  ONR/Pasadena 

1  ONR/San  Francisco 

2  NRL 

1  Code  2027 
1  Code  2629 

1  NORDA 

3  USNA 

1  Tech  Lib 
1  Nav  Sys  Eng  Dept 
1  B.  Johnson 

3  NAVPGSCOL 

1  Lib 

1  T.  Sarpkaya 
1  J.  Miller 

1  NOSC/Lib 

1  NCSC/712 

1  NCEL/131 

1  NSWC ,  White  Oak/Lib 

1  NSWC,  Dahlgren/Lib 

1  NUSC/Lib 


Copies 

11  NAVSEA 

1  SEA  05R24  (J.  Sejd) 

1  SEA  55W3  (E.  Comstock) 

1  SEA  55W33  (W.  Sandburg) 

1  SEA  55W31  (W.  Louis) 

1  SEA  55W31  (G.  Jones) 

1  SEA  55N2  (A.  Paladino) 

1  SEA  55N1  (S.G.  Wieczorek) 
1  SEA  63R31  (T.  Peirce) 

1  SEA  56X12  (C.R.  Crockett) 
1  SEA  62R41  (L.  Pasiuk) 

1  SEA  99612  (Library) 

1  NAVFAC/032C 


1 

NADC 

1 

NAVSHIPYD 

PTSMH/Lib 

1 

NAVSHIPYD 

PHILA/Lib 

1 

NAVSHIPYD 

NORVA /Lib 

1 

NAVSHIPYD 

CHASN/Lib 

1 

NAVSHIPYD 

LBEACH/Lib 

2 

NAVSHIPYD 

MARE 

1  Lib 

1  Code 

250 

1 

NAVSHIPYD 

PUGET /Lib 

1 

NAVSHIPYD 

PEARL/Code  202.32 

1 

NAVSEC,  NORVA/ 6660 . 03 ,  Blount 

12 

DTIC 

1 

AFOSR/NAM 

1 

AFFOL/FYS, 

,  J.  Olsen 

2 

MARAD 

1  Div  of  Ship  R&D 
1  Lib 


Preceding  Page  Blank 
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Copies 

1  NASA/HQ/Lib 

3  NASA/Ames  Res  Ctr,  Lib 

1  D .  Kwak 
1  J.L.  Steger 
1  Lib 

2  NASA/Langley  Res  Ctr 

1  Lib 

1  D.  Bushnell 
1  NBS/Lib 

1  LC/Sci  &  Tech 

1  DOT/Lib  TAD-491.1 

1  MMA 

1  National  Maritime  Res  Ctr 
1  Lib 

3  U  of  Cal/Dept  Naval  Arch, 

Berkeley 
1  Lib 

1  W.  Webster 
1  R.  Yeung 

2  U  of  Cal,  San  Diego 

1  A.T.  Ellis 
1  Scripps  Inst  Lib 

1  U  of  Cal,  Santa  Barbara/Tulin 

4  CIT 

1  Aero  Lib 
1  T.Y.  Wu 
1  A.J.  Acosta 
1  D.  Coles 

1  California  State  University/ 

Cebeci 

1  Catholic  U  of  Amer/Civil  & 

Mech  Eng 

l  Colorado  State  U/Eng  Res  Ctr 

1  Cornell  U/Shen 


Copies 

2  Harvard  U 

1  G.  Carrier 
1  Gordon  McKay  Lib 

1  U  of  Illinois/J.  Robertson 

4  U  of  lows 

1  Lib 

1  L.  Landweber 
1  V.C.  Patel 
1  C . J .  Chen 

1  Johns  Hopkins  U/Lib 

4  MIT 

1  Lib 

1  J.R.  Kerwin 
1  T.F.  Ogilvie 
1  J.N.  Newman 

2  U  of  Minn/St.  Anthony  Falls 

1  Lib 
1  Arndt 

1  U  of  Mich/NAME/Lib 

1  U  of  Notre  Dame/Eng  Lib 

1  New  York  U/Courant  Inst/Lib 

3  Penn  State 

1  B.R.  Parkin 

1  R.E.  Henderson 

1  ARL  Lib 

1  Princeton  U/Mellor 

1  U  of  Rhode  Island/F.M.  White 

1  Science  Application,  Inc. 

Annapolis,  MD 
C.  von  Kerczek 

1  SIT/Lib 

1  U  of  Texas/Arl  Lib 

1  Utah  State  U/Jeppson 
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Copies 

1 


Copies 

1 


Exxon  Math  &  System,  Inc. 


-  . 


Southwest  Res  Inst 
1  Applied  Mech  Rev 

3  Stanford  U 

1  Eng  Lib 

1  R.  Street,  Dept  Civil  Eng 
1  S.J.  Kline,  Dept  Mech  Eng 

1  Stanford  Res  Inst/Lib 

1  U  of  Virginia/Aero  Eng  Dept 

1  U  of  Washington/Arl  Tech  Lib 

2  VPI 

1  Dept  of  Mech  Eng 
1  J.  Schetz,  Dept  Aero  & 
Ocean  Eng 

2  Webb  Inst 

1  Lib 
1  Ward 

1  Woods  Hole/Ocean  Eng 

1  Worchester  Pi/Tech  Lib 

1  SNAME/Tech  Lib 

1  Bell  Aerospace 

1  Bethlehem  Steel/Sparrows  Point 

1  National  Science  Foundation/ 

Eng  Div  Lib 

1  Bethlehem  Steel/New  York/Lib 

4  Boeing  Corapany/Seattle 

1  Marine  System 
1  P.  Rubber t 
1  G.  Paynter 
1  H.  Yoshihara 

1  Bolt,  Beranek  &  Newman/Lib 

1  Cambridge  Acoustical 

Associates,  Inc. 

1  Exxon,  NY/Design  Div/ 

Tank  Dept 


1  General  Dynamics,  EB/ 

Boatwright 

1  Flow  Research 

1  Gibbs  &  Cox/Tech  Info 

2  Grumman  Aerospace  Corp 

1  Lib 

1  R.E.  Melnik 
1  Hydronautics/Lib 

1  Lockheed,  Sunnyvale/Waid 

1  Lockheed,  California/Lib 

1  Lockheed,  Georgia/Lib 

2  McDonnell  Douglas,  Long  Beach 

1  T.  Cebeci 
1  J.L.  Hess 

1  Newport  News  Shipbuilding/Lib 

1  Nielsen  Eng  &  Research 

1  Northrop  Corp/Aircraf t  Div 

1  Rand  Corp 

1  Rockwell  International 

1  B.  Ujihara 

1  Sperry  Rand/Tech  Lib 

1  Sun  Shipbuilding/Chief 

Naval  Arch 

1  TRW  Systems  Group/Lib 

1  TRACOR 

1  United  Technology/East 
Hartford,  Conn 

2  Westinghouse  Electric 

1  M.S.  Macovsky 
1  Gulino 
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CENTER  DISTRIBUTION 


Copies  Code  Name 

1  012.3  D.D.  Moran 

1  1500  W.B.  Morgan 

1  1504  V.J.  Monacella 

1  1506  S.  Hawkins 

1  1508  R.  Boswell 

1  152  W.C.  Lin 

1  1521  W.  Day 

1  1522  C.C.  Hsu 

1  1522  M.B.  Wilson 

1  1522  C.H.  Sung 

1  154  J.  McCarthy 

1  154  P.  Granville 

1  1540.1  B.  Yim 

1  1540.2  R.  Cumming 

30  1542  T.T.  Huang 

1  1542  N.C.  Groves 

1  1542  G.S.  Belt 

1  1542  Y.T.  Lee 

1  1542  M.S.  Chang 

1  1542  R.W.  Burke 

1  1544  F.  Peterson 

1  156  D.S.  Cieslowski 

1  1561  G.  Cox 

1  1564  J.  Feldman 

1  1606  T.C.  Tai 

1  1802.1  H.  Lugt 

1  1840  J.  Schot 

1  1843  H.  Haussling 

1  19  M.M.  Sevik 

1  1905.1  W.  Blake 

1  194  F.S.  Archibald 

10  5211.1  Reports  Distribution 

1  522.1  TIC  (C) 

1  522.2  TIC  (A) 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


1.  DTNSRDC  REPORTS,  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEM)  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 
BASIS. 


